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Abstract. A technique for detecting the presence of cloud in the field of view of an 
atmospheric Cerenkov telescope using a far infra red radiometer is described. Models for 
the radiative emission from clear and cloudy skies are tested and found to represent the 
measurements. 
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1. Introduction 

Measurement of atmospheric Cerenkov radiation is an important technique 
finding increasing application in high energy astrophysics (Fegan, 1997). 
Observations are customarily made under conditions of cloudless and moon- 
less dark skies. A long standing problem has been the need to monitor the 
clarity of the atmosphere — the radiating medium of the Cerenkov light 
production and detection system. A number of possibilities exist; recent 
suggestions, triggered by plans for future large air shower recorders, employ 
active systems involving laser probing of the atmosphere and measurements 
of the scattered light (Abu-Zayyad et al., 1997) or the monitoring of a com- 
plete hemisphere using a matrix of fixed far infra red (FIR) radiometers 
(Bird et al., 1997). Alternatively, monitoring of the brightness of a star in 
the field of a gamma ray telescope is a passive method already employed 
for ground based gamma ray astronomy (Armstrong et al., 1998). However, 
this method has not found general application because measurements of faint 
stars using inexpensive CCD detectors tend to be noisy giving an imprecise 
estimate of atmospheric clarity. Often the performance of the gamma ray 
telescope when measuring the Cerenkov radiation produced by the cosmic 
ray background is used as a monitor of the transparency of the atmosphere. 
This has the disadvantage of depending on the performance of the Cerenkov 
detector itself. 
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We report here a technique involving measuring the radiative temperature 
of the night sky and sensing the presence of cloud and other obscurants 
against the cold clear sky. The method is quick, sensitive and independent 
of the performance of the atmospheric Cerenkov detector. 

We present a range of observational radiometric data for the temperature 
of clear and cloudy skies and describe these data using two simple models 
for the temperature of clear sky and of clouds at various heights. 

2. Experimental Equipment 

Measurements of the sky temperature have been made with FIR radiometers 
(Heimann model KT 17 and KT 19) sensitive in the 8-14 jj,m waveband 
which coincides with an atmospheric window. The radiometers have a tem- 
perature range — 75°C — hlOO°C, a temperatue resolution of ±0.2°C, and 
an aperture of 2° defined by a germanium lens. The radiometers were cali- 
brated including allowance for the effects of the instrument temperature on 
the performance of the radiometer (~ 1°C change in indicated temperature 
for a 10°C change in the temperature of the instument). 

Measurements were made with radiometers mounted on the Mark 6 gam- 
ma ray telescope (Armstrong et al., 1998) in Narrabri, NSW, Australia and 
on steerable alt-azimuth platforms and on fixed mounts in Durham, Eng- 
land. Data from the measurements made at various angles to the zenith, 
throughout the day and night during different seasons over a period of years 
are available. Cloud base measurements were made in 1998 using a LIDAR 
ceilometer for a subset of the data. Full meteorological data are available for 
all observations. 



3. Results 

3.1. Correlation between telescope background count rate 
and sky radiation temperature 

The performance of a gamma ray telescope and the stability of the atmo- 
sphere has often been monitored on the basis of the background cosmic ray 
count rate of the telescope. This method has a number of disadvantages. It 
depends upon the integrity of the telescope performance and it is governed 
by the statistics of counting (at rates with current telescopes up to 1000 
cpm). We show in Figure |] the correlation between the background count 
rate of the telescope and the radiative temperature of the sky during an 
observation using the Mark 6 telescope at Narrabri. The sky temperature is 
seen to be very sensitive to the presence of the warmer clouds and correlates 
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Figure 1. The correlation between the background counting rate of the Mark 6 gamma 
ray telescope (solid line) and the radiative temperature of the sky (broken line). 



with the background count rate of the telescope. The presence of clouds 
caused a decrease in telescope count rate and an increase in sky tempera- 
ture. It thus presents an instantaneous, sensitive and independent monitor 
of the clarity of the sky in the beam of the gamma ray telescope. 

3.2. Zenith angle dependence of temperature of a clear sky 

We show in Figure || the measured temperature of a clear sky using the 
Narrabri radiometer as a function of zenith angle of observation in the range 
10° — 70°. The plot is based on data from observations during two nights when 
there was a small difference in air temperature, which is noted. The radiative 
temperature of the clear sky increases as the zenith angle is increased and 
the measurement is made through an increasing slant thickness of warm 
lower atmosphere. The data of the figure also show that the temperature of 
the clear sky from night to night is reproducible. 

3.3. Diurnal variation of temperature of a clear sky 

The profile of the sky temperature (measured at the zenith near Durham, 
England) during a typical 24 hour period during which the sky was clear is 
shown in Figure ||. The air temperature at screen height is also plotted and 
it is noted that the vertical clear sky temperature tracks the air temperature 
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Figure 2. The variation in the clear sky radiometric temperature as the zenith angle at 
which the measurement is made is increased from 0° to 67°. The data were obtained on 
two nights and there is a small variation in the air temperature between the observations 
at zenith angles < 30° and > 30°, as shown in the figure. 

at screen height. Such data taken over a long period are the basis for a simple 
phenomenological model linking the observed sky radiometric temperature 
with the air temperature and humidity at screen height — see Section [|. 

3.4. Annual variation of temperature of a clear sky 

The temperature of a clear sky measured at the zenith at 1400 UTC near 
Durham, England averaged over a four week period during July, September 
and November is shown in Figure |j. A systematic reduction in the sky 
temperature occurs during the transition from summer to winter illustrating 
the annual variation. 

3.5. Temperature of clouds 

Figure [B| shows measurements of the sky temperature at the zenith and 60° 
to the zenith during a 24 hour period when clear skies were replaced by total 
overcast clouds with a base at 600 m at about 1200 UTC. The temperature 
of the clear sky is seen to be about 40° C cooler than than that of the low 
overcast cloud which is very similar to the ground level air temperature. 
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Figure 3. The variation in the clear sky radiometric temperature (solid line) over a 24 hr 
period. The profile of the air temperature at screen height (broken line) is also shown. 




Figure 4- The variation in the clear sky radiometric temperature at the zenith at 1400 
UTC averaged over 4 week periods between from July and November. The data were 
recorded near Durham, England. 
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Figure 5. The sky radiative temperatures at the zenith (broken line) and at 60° to the 
zenith (solid line) in intervals of cloud and clear sky during a 24 hour period. 

3.6. Variation of radiometric temperature of clouds with 

HEIGHT 

The temperatures of clouds with bases at different heights were measured 
at the zenith using a radiometer boresighted with a LIDAR ceilometer - 
see Figure ^. The radiometer has an aperture of about 2° and on occasion 
a part of the cold clear sky is viewed when the cloud base indicated by the 
narrow angle LIDAR is low. This is the origin of some spread of the data in 
Figure |§ An increase in cloud base height from 1000 to 7000 ft corresponds 
to a decrease in measured cloud temperature from +10°C to 0°C. 

4. A Model for the Temperature of Clear Skies 

4.1. A SIMPLE EMPIRICAL MODEL 

The data of Figure [7| show the relation between the radiometric temperature 
of a clear sky and the air temperature at screen height for measurements 
near Durham, England. They suggest that the major factor governing the 
clear sky temperature is the temperature of the air at the observation level. 
This is representative of the temperature of the atmosphere through which 
the measurement is made. 
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Figure 6. The measured cloud radiative temperature plotted against cloud base height 
The contours specify the frequency of observations of a given cloud radiative temperat ure 



and base height combination. The line represents the model prediction — see Section 5.2 




Figure 7. The relationship between the radiative temperature for a clear sky observed at 
the zenith and the air temperature at screen height. 
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An approximate representation of the sky radiative temperature is given 

by 

T sky = 1.4T - 66.2 (1) 

where T s k y is the sky radiative temperature (in °C) and To is the air tem- 
perature at screen height (in °C). 

4.2. The model of Idso (1981) 

A number of models for the emission of long wave radiation from clears skies 
have been developed, normally for use with the modelling of meteorological 
conditions and in energy-balance models (see e.g. Monteith and Unsworth, 
1990). Because of this principle application, the majority of these models 
predict the longwave radiation integrated over a complete hemisphere and 
not from a small area of the sky. A recent example of this type of model is 
that due to Prata (1996). 

Idso (1981) described a model which can be used to predict the radiation 
from a small area of clear sky. The model incorporates both the effects of 
air temperature and vapour pressure at ground level in the prediction of sky 
temperatures. 

Idso showed that the full spectrum radiated by a clear sky, i2 a , may be 
represented by 

R, = e a aT$ (2) 

where e a is the full spectrum clear sky emittance, To is the screen level air 
temperature in K and a is the Stefan-Boltzmann constant. 

An equation for e a , appropriate to the 8-14 /im wave band at the zenith, 
was given by Idso as 

e a = 0.24 + 2.98 x 10 _8 eg exp(3000/T ) (3) 

where eo is the screen level (ground level) vapour pressure and can be cal- 
culated using an equation given by Wiesner (1970). 

4.3. Validation of the Idso model 

Data were selected from observations on clear days and the vapour pressure 
and air temperatures measured on these days were employed as input to 
the model to calculate the radiometer sky temperature. The predicted sky 
temperature was compared with the measured values. Figure || shows the 
measured sky temperatures and the predictions of the model for varying 
air temperatures at screen height (the variation in vapour pressures were 
also allowed for). The predicted and observed variation of temperature can 
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Figure 8. A comparison of the predictions of the Idso model (•) with experimental data 
(o) for the variation of radiometer temperature for clear skies with air temperature at 
screen height. The solid line is a linear fit to the experimental data and the broken line a 
linear fit to the modelled points. 

each be fitted by a straight line. It can be seen that the predicted and 
observed temperatures represented by the straight lines differ by 2°C at low 
air temperatures and 10°C at high air temperatures. 

This simple model therefore provides a reasonable representation of tem- 
perature of the clear sky. Explanations may exist for the differences noted. 
For example, it is possible that the constants of the model should be adjust- 
ed to take account of the difference in climate in the northeast of England 
(where our measurements were made) and in Arizona (where the original 
data used to validate the Idso model were recorded). A small adjustment 
has been made to equation [3| specifying the 8-14 /j,m clear sky emissivity 
(e a ) which reduces the values by ~ 10% and reconciles the predictions of the 
model and our measurements. The original formulation by Idso has been 
modified to 

e a = 0.22 + 2.98 x 10~ 8 eg exp(3150/T ). (4) 

This is the form of the equation appropriate for skies observed under 
climatic conditions found in northeastern England. Figure || shows the com- 
parison of the predictions from this modified model and experimental mea- 
surements; the fitted lines are indistinguishable. 
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Figure 9. A comparison of the predictions of the modified Idso model (■) with exper- 
imental data (□) for the variation of radiometer temperature for clear skies with air 
temperature at screen height. The predictions (broken line) and the observations (solid 
line) are now indistinguishable. 



5. A Model for the Temperature of Clouds 

If clouds are present in the sky, they produce an additional contribution to 
the thermal radiation detected at ground level. For a cloudy sky, therefore, 
the measured sky radiative temperature is increased. As is the case for clear 
skies, much of the discussion in the literature has been of models for the 
prediction of the thermal emission from a complete hemisphere for use in 
meteorological models (e.g. Alados-Arboledas et al., 1995; Malek, 1997). 
However Kimball et al. (1982) have developed a model that can be used to 
predict the cloud temperature at a specified zenith angle. 

5.1. The model of Kimball et al. (1982) 

This model for the emission from a cloudy sky is based on the clear sky model 
of Idso (1981) but with an additional term to include radiation effects of the 
cloud. In this model the contribution to the thermal radiation from clouds, 
R c , is 

R c = rAeJaT* (5) 

where e c is the cloud emittance (a value dependant on cloud type in the 
range 0.5 to 1.0), / is the fraction of black body radiation emitted in the 8 
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- 14 /j>m band, (taken as 1.0 in this work), A is the fraction of sky covered 
with cloud, a is the Stephan-Boltzmann constant and r is the transmittance 
of the atmosphere in the 8-14 fim window where 

t = 1 - e z (6) 

and e z is the 8-14 /im sky emittance (Idso, 1981) — see equation |3[ 

T c is the cloud temperature, which is calculated assuming the standard 
temperature lapse rate (see e.g. Lutgens and Tarbuck, 1998). According to 
this the equation for the temperature of the cloud is 

T c = T - 0.0065Z (7) 

where To is the air temperature at screen height in K and Z is the height of 
the cloud in meters. 

The total predicted sky thermal radiation, R^, received at the earth's 
surface from the cloudy sky is then 

R T = R a + R c (8) 

where R a is the clear sky radiation — see equation |2[ 

5.2. Validation of the model 

The cloud model developed by Kimball et al. (1982) was used to predict 
radiometer temperatures for cloudy skies which may be compared with 



observations. Figure |10| shows a comparison between the predictions of the 
model and observations on a day which included changes in cloud cover. 
The data were taken in August in northeast England. The best fit of the 
model predictions and observations (which is shown) was when the cloud 
emissivity was taken to be 0.8 (a reasonable mid-range value). 

Comparison of the predictions of the Kimball et al. model with the mea- 
surements of temperature for clouds at different heights are shown in Figure 
§. There is good agreement between the solid line representing the predic- 
tions and the contours repesenting the data when a value of 0.8 is used for 
the cloud emissivity. 



6. Summary 

The use of a FIR radiometer sensitive to the 8-14 /im waveband mounted 
on a ground based gamma ray telescope provides a simple and inexpensive 
method of reliably monitoring the atmospheric clarity in the telescope's field 
of view. The response of the radiometer is prompt and is independent of the 
performance of the telescope itself. 
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Figure 10. A comparison of the predictions of the Kimball et al. model (solid line) with 
observed cloud/sky temperature (broken line) during a changeable day. 



The predictions of simple semi-empirical models of the longwave emission 
from clear skies (Idso, 1981) and cloudy skies (Kimball et al., 1982) have been 
compared with our measurements. They provide theoretical underpinning 
for the observed variation of radiative temperature of the sky and clouds 
with the air temperature and relative humidity at screen height. 
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